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Abstract. An evolution scheme is developed, based on Sorkin 
algorithm, to evolve the most complex regular tridimensional polytope, 
the 600-cell. The solution of 600-cell, already studied before, is 
generalized by allowing a larger number of free variables. The 
singularities of Robertson- Walker (RW) metric are studied and a reason 
is given why the evolution of the 600-cell stops when its volume is still 
far from zero. A fit of 600-ceirs evolution with a continues metric is 
studied by writing a metric generalizing Priedmann's and including the 
600-cell evolution too. The result is that the 600-cell meets a causality- 
breaking point of space-time. We also shortly discuss the way matter is 
introduced in Regge calculus. 
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1. Introduction 

Regge calculus is 40 years old and has been extensively studied. In its original form it 
applies only to an empty spacetime, and the problem of how to introduce matter has 
always been a major one. 

A key test of Regge Calculus in presence of matter which has been studied 
many times in literature is the Friedmann universe of dust. The first papers on 
the Friedmann universe did not use simplicial complexes. The authors made use 
of non simplicial blocks for which they had to give some other information besides 
the lengths of the edges. This information was deduced by symmetry. Collins and 
Williams (1973) used three tridimensional regular polytopes to describe the spatial 
sections with topology of the Friedmann universe of dust. These authors assumed 
g, the dust density, to be constant between two spatial slices. Brewin (1987) wrote in 
the discrete formalism of Regge Calculus the expression of the action. His expression 
is equivalent to the one that would result by assuming a continuous density ^ oc 

Besides the regular polytopes, Brewin also used other non regular tridimensional 
polytopes. The evolution of all his polytopes, using Regge equations, stopped before 
the volume of spatial 3-section became zero: a point was reached where the equations 
had no solution. He saw that this point approaches the limit value (spatial section 
of null radius) as the polytope approaches the 3-sphere of the continuum solution. 
To overcome the obstacle, Brewin continued the evolution reverting to differential 
equations. 

The first authors who used a 4-dimensional simplicial complex for the Friedmann 
universe of dust were Barrett et al (1997). They evolved a spatial section made of a 
regular 3-polytope having n vertices, and "discretized" matter by placing one mass 
point on each timelike edge joining a vertex i of the polytope with its evolute i' of the 
next spatial section. Thus they found the action: 



where the first sum is over the bones h of the manifold and Si, is the defect of the bone 
h with area Ai,. The second sum is done over every vertex i of the first spatial section. 
For a regular polytope having n = 120, a mass M/120 is put in each vertex i. 
The Regge equations are found by extremizing the action (1): 



where the sum is done over all vertices k joined to the edge [ij] of length 1^. The bone 
[ijk] has Sijk as defect and A^j^ as area. They made use of the evolutive scheme found 
by Sorkin (1975) by assuming that the vertices of the 600-cell could be subdivided in 
four classes each composed by vertices not joined together. 

We have found that the vertices of the 600-cell actually belong to five classes and 
have investigated if this difference might lead to different results. 
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All evolutive schemes described above exhibit a common behaviour: before 
reaching the "true" singularity of Friedmann universe, the evolution stops because 
Regge equations have no solution. It is not clear if this happens because of the 
discretization intrinsic in Rcggc calculus, or for other reasons (see section 5). We 
decided to investigate this point in greater detail. 

This paper is organized as follows. In section 2 we sketch the results of our 
calculations. In section 3 we embed the closed Robertson-Walker metric into a five- 
dimensional Lorentzian space. We find it useful in order to understand the problems 
related to the evolution of our lattice. The nature of the singularities belonging to 
the closed Robertson-Walker metric is discussed in section 4. In section 5 we propose 
a deeper reason why our simplicial approximation of the Friedmann universe ceases 
to exist before the volume of the spatial 3-section vanishes. By fitting the simplicial 
solution with a continuous metric (which has the same behaviour with respect to 
the reaching of singularities) we find that the stop condition is caused by a novel 
singularity of the metric (section 6). Finally, our conclusions are reported in section 7. 

2. Results of the evolution 

The vertices of the 600-cell belong to five classes each including only vertices not joined 
together. Using the classification given by Coxeter (1973), one of these five classes is 
composed by the following vertices: 

— {Aqj Aio, A20, ^407 ^507 Bn, B27, -B37, -B47, -B57, 
Ci, Cii, C21, C31, C41, C51, Dg, Dig, D28, -D38, D48, D^g}. 

The other classes (3, 7, 5, e are obtained by adding 2, 4, 6, 8 to the subscripts. 

It should be noted that every evolution consists in solving 5 systems of 4 equations 
in 4 unknowns each, given the lapse and shift freedom. Furthermore the simplicial 
sandwich does not have the same grade of symmetry as the initial spatial section, 
where all vertices were equivalent. So the Sorkin evolutive algorithm we used breaks 
the symmetry among the vertices a, (3, 7, 5, e. 

We have generalized the evolution of the 600-cell by imposing that the edges of 
the same type, like [a/?], were all equal. In this way only four unknowns and four 
equations remained when each vertex was evolved. By allowing the edges joining 
vertices belonging to two different classes to be possibly different one from another, 
we introduce new degrees of freedom in the lattice. We have verified that this 
enhancement of freedom causes no problem in the evolution. 

We performed the evolution of a 600-cell starting from the instant of time 
symmetry. Our results are in figure 1. The evolution stops at a finite value R^i 
of the radius R (the radius of the three-sphere equivalent to the 600-cell) which is 
approximately 1/4 of the initial value (i?o = 4.244...). For R < i^m the Regge 
equations give no solution. Further numerical details may be found in De Felice and 
Fabri (2000). 
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3. The embedding 5-space 



We shall find convenient to embed a space-time with closed Robertson- Walker metric 
into a 5-dimensional Lorentz space. This will be of help for understanding the 
behaviour of the Regge equations when applied to the evolution of a 600-cell space 
section. 

The metric under consideration can be written in this way: 

ds^ = -dt^ + [dx^ + sin^ x (di?^ + sin^ ^ d(^^) ] , (2) 

where R depends only on the universal time t. 

The manifold with metric (2) can be embedded into a 5-dimensional space endowed 
with the (1+4) metric 

da^ = -dv^ + dw'^ + dx'^ + dy'^ + dz'^. (3) 

In the following we will refer to v as "outer time." 

The metric (3) can be rewritten using coordinates v, R, x, (f by the substitution 

w = R cos X 

X = R sin X sin cos (p 

(4) 

y = R sin X sin -d sin (p 
z = R sin X cos 

Then 

da^ = -dv^ + dR^ + R^ [dx' + sin^ x (di$>' + sin^ ^dif'')]. (5) 

A 4-dimensional submanifold of (5) can be assigned by imposing that i? is a function 
only of the outer time v: then we have {R = dR/dv) 

ds^ = - (1 - R^) dv^ + R^ [dx^ + sin^ x (d^^ + sin^ ^ dip^)] . (6) 

In order that equations (2) and (6) may coincide it is necessary that 

dt" = (1 - R?) dv\ (7) 

from which it follows that R^ < 1. We are interested in the case where R{v) is an 
even function, so that v = is an instant of time symmetry. Integrating equation (7), 
t can be found as a function of v. 
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4. Singularities 



In order to study the singularities of metric (6) let us compute the quadratic 
Riemann invariant: 



12 



1 

+ ^ 



(1 - Wf 



(8) 



We will assume as a necessary and sufficient condition for a singularity that Q becomes 
infinite. A glance to equation (8) shows that the singular points are those values of v 
where \R\ — 1 ori2 = 0or \R\ = +oo. Let us discuss the possible cases. 

1. The points R = 0. These coincide with the infinite contraction of the universe 
when the volume of the spatial section goes to 0. In these points the density 
reaches infinite values. By a volume-vanishing (VV) singularity we will mean 
one of these points. 

2. The points R = ±1. If we are not in the first case we can say that the radius 
of the space section increases or decreases with speed tending to that of light. If 
one tries to extend the manifold beyond these points, and if increases further, 
the metric becomes positive defined. So it can no longer represent a manifold of 
General Relativity. This is tantamount to saying that the spacelike sections are 
not causally connected: the points of two different spatial sections can be joined 
only by spacelike geodesies. These points will be said to be causality-breaking 
(CB) singularities. Actually, in a rigorous sense a CB singularity is a border of 
the manifold, and the extension is meaningless. 

3. The points R = ±oo. If we are not in the previous two cases the singular 
behaviour is shown by the geodesies deviation, which becomes infinite. We will 
refer to these points as geodesies-deviation (GD) singularities. 

It should be noted that in general the GD points do not represent a singularity 
according to the definition used in the context of Hawking's theorem, see for example 
Wald (1984), where the presence of a singularity means that the universe is geodesically 
incomplete. An example of a GD point which does not stop any geodesic is the point 

V — for the universe having R[v) — Rq + Av*^^. 

These singularities can also occur together in the same point: an example is the 
hypothetical universe whose radius has the form R{v) = v — Av^l"^ (the critical point 

V = is a VV-CB-GD singularity). Another interesting case of singularity merging is 
the Friedmann metric discussed in subsection 4.2 below. 



Jf..l. Matter-limited singularities 

From Einstein equations for the closed Robertson-Walker metric we have for the 
diagonal components of T^js in an orthonormal basis: 



6 



A De Felice and E Fabri 



3 

Q = 



(9) 

1 2RR+1-R^ 

It is well known that the stress-energy tensor must satisfy 
Q > \p\- 

In our case 

_2<-5^<l. (10) 

In general it can happen that condition (10) be not satisfied for values of v for which 
the metric still has a "geometrical sense." In fact a CB or GD or CB-GD singularity 
does not satisfy equation (10). Instead a VV point presents no problem. It should 
be noted that we cannot become aware of these new "physical" matter-limited (ML) 
points, wherein g = \p\, by looking only at the expression of the metric. 

4-2. Friedmann metric 

Since dt = Rdrj, from equation (7) we obtain 



v = J [R^ -{dR/driff'dr]. 



For the Friedmann metric we have: 

77 

V = 2i?o sin - 

and 

2 

R = Ro cos^ - = i?o ^ 



2 4i?o 
Metric (6) takes the following form: 

ds^ = -^ dv^ + R^ [dx^ + sin2 X {^^^ + sin^ d^^)] . 
Rq 

Calling Vc the value of v where R = —1 and the positive value where 
i? = 0, for the Friedmann metric these two singularities occur in the same point, 
i.e. — Vc — 2Rq. So, in the 5-space, in correspondence of the big crunch 
the speed of contraction equals the speed of light, i.e. the big crunch is a VV-CB 
point. This behaviour is shared by all matter satisfying p = w g, where is a 
constant, — 1/3 < < 1, since then R^ = 1 — R^'^^'^, A a constant. Seeing 
that R = ~A^{l + 3w)R^'^,m Zel'dovich's interval, < < 1, no GD singularity 
is present besides the other two. 



Singularities of RW metric 



7 



5. The problem of the 600-cell 

What is the nature of the critical point for the 600-cell model? Does it reach a 
singularity, and in this case which of them? Brewin (1987) gives a partial answer to 
these questions by noting that the following relation occurs: 

Al ^ 

1, 

T 

where r is the interval of universal time between two consecutive spatial sections 
and Al is the difference of the lengths of their edges. The collapse becomes so fast 
that the vertical edges are to become spacelike. Actually in our model the above 
fraction exceeds 1 and reaches a not easily interpretable value (approximately 2.6). 
Barrett et al (1997) saw the position of the stop point is quite independent of the 
timelike interval between two consecutive spatial sections and advanced the idea that 
the evolution should be continued by making spacelike the edges of the lattice which 
were timelike before the stop. However they did not pursue the idea. So what is 
the true meaning of the above condition, and how can we continue the iterations (if 
possible)? 

In order to answer all these questions let us now reconsider equation (7). For each 
evolutive step we know At (= r) and AR. Then we can find an approximate value of 

1 /9 

Av by calculating Av — (At^ + Ai?^) and compute as a function of v (figure 2). 
We also give AR/ Av dR/dv as a function of R in figure 3. We can see that this 
derivative in the last iterations tends, in absolute value, to the unity. In this sense 
the 600-cell meets a CB singularity before its volume vanishes. So we expect that for 
the 600-cell Vc < Vq. This leads us to believe that the evolution cannot be continued 
beyond this point. 

We can get a clearer understanding of the situation examining a case easier to be 
studied in detail. Consider the following partial differential equation: 



^ + ^ = -t-xt-lx^, with y(x,0) = i. (11) 

Its solution is 



y = l-lt'-'^tx\ (12) 

We can see from equation (12) that the relation y = is satisfied along two curves. 
There are also two curves (actually two parabolas) where dy/dt = ±1. These four 
curves are tangent two by two in two points, x — and t — ±1, wherein both conditions 
hold. 

Let us try to solving equation (11) by using a numerical approximation. If we use 
Lax method (see Press et al (1992)), taking h as the spatial a;-step and k the t-step 
(h > k for the Courant condition), the above differential equation becomes 

Vm+ljU ym—l,n '^ym,n+l (i/m+l,n ~l~ ?/m— l,n) , j. 1 2 

2h 2k ~ ~ ~ " ~ 2 ^m- 
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The solution can be easily worked out and we can see that the new four curves are 
not tangent any longer. Furthermore the region where \dy/dt\ < 1 is strictly included 
into the region y > 0. So the two points in a; = are splitted and for every x the 
solution starts and ends in two "CB points" {dy/dt = ±1). This discrepancy vanishes 
as h, the x-step, reaches smaller and smaller values. 

The above argument shows that our 600-cell model does not stop because it is going 
to reach the singularity R = — a situation which, like in all numerical methods, could 
give rise to problems. Instead our model stops because it fits a universe that reaches 
another kind of singularity {R— —1) 



6. Study of a "new" class of solutions of Einstein equations 

We have just seen that the behaviour of the evolution for the 600-cell looks somewhat 
different from the Friedmann solution, in its reaching a CB singularity before the 
vanishing of R. So it appears expedient to study a generalisation of Friedmann metric, 
also having that property. For the 600-cell a good empirical fit of ) is given by 

R{v) = Ro- with w 1.128, 

ARo 

whereas a = 1 is the Friedmann metric. 
It is useful to re-define 

def 2i?o 



def 2i?o 

a 



Vn = = avr. 



so 



Vc 

• For < a < 1 we have that —v^ < v < Vq and the singularity is of VV-type. 

• For a > 1 we have that —Vc < v < Vc and the singularity is of CB-type. This is 
the case our model fits. 

• For a = 1 the solution is the Friedmann metric. So this case is a watershed 
between two different behaviours of the general metric. 

It is not difficult to find the expression of t as a function of v: 

1 /2 

t = I Vcarcsin(v/vc) + f'*^ [l — (v/vcY] 

and to show that the graph of R{t) is still a cycloid, like for Friedmann metric, but 
scaled in i? by a factor 1/a^ and translated upwards by | (a^ — 1). 
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If a < 1 this condition causes no problem because v"^ > Vq. If a > 1 the ML point v 
is the first one the metric meets in its evolution. In fact 



Note that in order to have -D^ > the value of must not exceed 4. So the conditions 
for the stress-energy tensor reduce the interval of definition ofato0<a<2. Ifa<l, 
the metric ends in a VV singularity reached the more softly the smaller the value of 
a. If a > 1, the metric reaches a ML point before the CB one. The Friedmann metric 
(a = 1) is more "pathological" because its VV singularity is also a CB point, since 
% = Vc- 

For this metric it is also possible in principle to write down an equation of state 
of matter 



7. Conclusions 

We have seen that the evolution of the 600-cell does not describe the Friedmann 
universe well. Instead we can think of it as the evolution of a different type of matter. 
But where does such a difference come from? In our opinion, one should consider the 
following three points as playing an important role in answering this question. 

First, the 600-ccll is a coarse approximation to a 3-sphere. For every approximating 
method, Regge Calculus included, the smaller the step interval the better the fit. When 
this interval reduces, the two solutions (numerical and analytical) tend to coincide. 
In Regge Calculus, in order to get a better approximation, one should increase the 
number of tetrahedra of each space section, i.e. the 600-cell should be substituted by 
another non-regular polytope nearer to a 3-sphere (Brewin (1987)). 

The second aspect is deeper since it is closely related to the grounds of the 
Regge Calculus itself. We have seen that Regge equations give rise to a different 
evolution of the universe, ending in a CB point. This happens because they are only 
an approximation of Einstein equations. 

Let us discuss this point by analogy with a related topic. In solving differential 
equations it is usual that a numerical method works better than another. For example 
let us take the following equation 



~2 2 12 



{a' - 1) < vl 




but of course we do not attach any physical meaning to this result. 



da; 
dt 



f{x,t) 
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and try to find a solution. If we used Euler's method we should write 
whereas, using the central-differences method, we would obtain 



Xh+i — Xk-i dx , , d^x , , dx 



So, even if both methods lead to the correct solution when h approaches to 0, the 
central-differences method is generally more accurate and converges more quickly than 
Euler's. According to the approximating method we choose, we have a smaller or 
greater truncation error. Furthermore, given that an approximate value for dxjdt 
implies the presence of other additional derivative terms of order greater than one, 
solving a differential equation through a numerical method is more like to solve another 
differential equation, which, as we have seen in section 5, of course has different 
solutions which may behave qualitatively different from the one we are looking for: 
for instance, the approximate solution may exhibit singularities not belonging to the 
correct solution. 

The third point is related to the way matter is taken into account. We followed 
Barrett et al (1997) by putting each mass point along the vertical edges, i.e. the "rest" 
geodesies that pass through the vertices of the simplicial complex. This appears to 
be the only easy way to introduce matter in a simplicial complex; so doing, however, 
a further discretization is introduced, not logically related to the original Regge idea. 
Alternatively, the mass points could be placed along other "rest'' geodesies; then the 
proper lengths would be a function not only of r but also of other edges belonging to the 
region between two spatial sections. In this case we should obtain a novel expression 
for the action and of course new Regge equations which might behave differently. 
Furthermore, we could put more than one single mass point in each tetrahedron, or 
even think of tetrahedra uniformly filled with dust. It is not clear how this change 
could reflect on the behaviour of the model. 

The undesired halting of the iteration scheme appears to be a property of all 
approaches attempted thus far to the Friedmann universe of dust, irrespective of their 
different treatment of matter. It is then reasonable to ascribe such behaviour either to 
the spatial discretization or to the Regge equations themselves, as explained before. 
Our opinion — admittedly not entirely justifiable — inclines towards the latter. 
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1 2 3 4 5 6 7 



t (universal time) 

Figure 1. Evolution of the radius i? as a function of the universal time t through 
Regge Calculus and General Relativity 
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Figure 2. The behaviour of i? as a function of v. 




Figure 3. The behaviour of |Ai?/Ai;| « \R\ as a function of v for the 600-cell 



